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1 THEORETICAL AND EXPERIMENTAL INVESTIGATION OF THE 


;ONIC-FLOW FIELDS BENEATH SWEPT AND UNSWEPT WINGS WITH 
hos TABLES OF VORTEX-INDUCED VELOCITIES 


By William J. Alford, Jr. 
SUMMARY 


The flow-field characteristics beneath swept and unswept wings as i! 
determined by potential-flow theory are compared with the experimentally . 
determined flow fields beneath swept and unswept wing-fuselage combina- ‘ 


cle lr wear et eee or: 
a 


7 Shas eri SST ee 


x . tions. The potential-flow.theory utilized considered both spanwise and : | 
: chordwise distributions of vorticity as well as the wing-thickness effects. : a 
ee The perturbation velocities induced by a unit horseshoe vortex are qnelused 

in tabular form. f + 


The results indicated that significant chordwise flow gradients 
existed beneath both swept and unswept wings at zero lift and throughout 
the lift range. The theoretical predictions of the flow-field character- 

= ‘istics were qualitatively correct in all cases considered, although there 
were indications that the magnitudes of the downwash angles tended to be ay 
overpredicted as the tip of the swept wing was approached and that the 
sidewash angles ahead of the unswept wing were underpredicted. The cal- a 
culated effects of compressibility indicated that significant increases 
in the chordwise variation of flow angles and dynamic-pressure ratios 
should be expected in going from low to high subsonic speeds. 
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INTRODUCTION 


\ The almost universal present-day employment of external stores, such 
i as missiles, bombs, or fuel tanks on fighter airplanes, and nacelles on 
Pee bomber airplanes, has indicated the need for more detailed information 
a regarding the flow characteristics in the vicinity of the wing in order 
fae to estimate the aerodynamic loads on these objects when fixed in the wing 
'¢. flow field and to evaluate the launching and jettison characteristics of 
missiles, bombs, or fuel tanks. In addition, numerous present-day air- 
od planes are incorporating wing sweep, lower aspect ratios, and shorter 
tail length, all of which may tend to bring the various airplane compo- 
nents in closer proximity to the wing. 
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, For airplane designs of the past, in which the component parts (for 
example, the wing and the tail) were separated by reasonable distances, i 

H the wing-interference effects could be calculated with sufficient accuracy " 

by a number of horseshoe vortices distributed along a single lifting line 

‘s (refs. 1 to h). However, because of the mathematically singular nature of 

. the single vortex, this ‘theory is valid only for regions that are at a 

1 distance of at least one wing chord from the vortex location. (See ref. 1.) 


The purpose of the present paper is to show that the flow character- 
istics beneath the wing can be calculated if the lifting wing is assumed 
to be represented by a mltiple arrangement (both chordwise and spanwise) 
of horseshoe vortices and if the effects of thickness are accounted for. 
The velocities induced by the airfoil-section thickness distribution, 
which are often neglected, are considered by using the appropriate sin- 
gularity (source sink) distribution (ref. 5) in conjunction with simple 
sweep theory (ref. 6). Detailed experimental flow fields were obtained 
around swept and unswept wing-fuselage combinations and are compared with 
the wing-~alone theoretical flow fields. 
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The details of the calculative procedure are developed in appendixes. 
The velocities induced by a unit horseshoe vortex in the chordwise, ver- 
tical, and lateral directions for a large range of distances are included 
in tabular form. The calculated first-order effects of compressibility 
on the flow characteristics for a subcritical Mach number of 0.80 are 
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also presented. 
: a 
[ SYMBOLS 
‘ A aspect ratio 
b wing span, ft 
c local wing chord, ft 
e mean aerodynamic chord, ft 
Cav average wing chord, ft 
: Cy wing-section lift coefficient 
— C1 section lift-curve slope 
‘- Cy, total lift coefficient ® 
CL incompressible lift-curve slope ‘. 
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i Cho, compressible lift-curve slope 


i Cp drag coefficient 

Sar 

~ Cn pitching-moment coefficient measured about 

mean aerodynamic chord 

1 fuselage length, 7.61 ft 

Me 5 wing area, sq ft i 

i s semiwidth of horseshoe vortex, ft : i 
Omax maximum fuselage diameter, 0.70 ft : ; 
me t airfoil thickness, ft 4: 
A taper ratio ‘ 
~ A local sweep angle, deg ' | 
7 Vv free-stream velocity, ft/sec : ' 


{ 
| VR resultant velocity, ft/sec 
K 


u backwash perturbation velocity in direction of x-axis, positive 
: rearward (fig. 3), ft/sec 


ec le wrt ror- re ocr 


Hl Ug backwash perturbation velocity induced by two-dimensional 
7 airfoil-section thickness distribution (see appendix A), 
f ft/sec 
Vv sidewash perturbation velocity in direction of y-axis, positive 
to the right (fig. 3), ft/sec 
w downwash perturbation velocity in direction of z-axis, positive 
downward (fig. 3), ft/sec 
a1 local dynamic pressure, 1b/sq ft 
| 
; Io free-stream dynamic pressure, 1b/sq ft 
€ downwash angle between free-stream-velocity vector and resultant- 


velocity vector in xz-plane, positive downward (fig. 3), deg 


0 sidewash angle between free-stream-velocity vector and resultant- 
velocity vector in xy-plane, positive toward left wing tip 
(fig. 3), deg 
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right-hand Cartesian coordinate system in which x is positive 
downstream, y is positive to the right, and z is positive 


upward (fig. 3), ft 


in the x-, y-, and z-directions, respectively, from 


s 
space point of interest to centroidal location of mth, nth 


vortex 


spanwise vortex index (see appendix A) . 

chordwise vortex index (see appendix A) 

inclination of wing from zero-lift attitude, deg 
three-dimensional vortex circulation strength, ft/sec 
two-dimensional vortex circulation strength, £t2/sec 
perturbation velocity potential, ft2/sec 


two-dimensional perturbation velocity potential (also referred 
to as chordwise accumulation of vorticity when increased by 


a factor of 2.0), ft@/sec 


backwash factor (see appendix B) 
sidewash factor (see appendix B) 
downwash factor (see appendix B) 


Mach number 


p =\/1 - Mw? 


Subseripts: 


a 


n 


additional or lift-indueed characteristics 


characteristics of airfoil section normal to local lines of 
constant percent thickness 


characteristics of streamwise airfoil section in two-dimensional 


flow 
characteristics referred to half-chord line 


characteristics referred to quarter-chord line 


characteristics referred to trailing edge 


Primes indicate equivalent incompressible characteristics. 
Bars indicate centroidal locations of the vortices. 
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MODELS AND TESTS 


The models about which the flow surveys were made consisted of both 
swept- and unswept-wing—-fuselage combinations. Drawings of the wing- 
fuselage combination are presented in figure 1. The wing of the swept- 
wing—fuselage combination had 45° sweep of the quarter-chord line, an 
aspect ratio of 4.0, a taper ratio of 0.3, and NACA 65A006 airfoil sec- 
tions parallel to the plane of symmetry. The wing of the unswept-wing—. 
fuselage combination had 0° sweep of the one-half-chord line, an aspect 
ratio of 3.0, a taper ratio of 0.5, and NACA 65A004 airfoil sections 
parallel to the plane of symmetry. The fuselage consisted of an ogival 
nose section, a cylindrical center section, and a truncated tail cone. 
The fuselage ordinates are presented in table I. 


The tests were made in the Langley 300 MPH 7- by 10-foot tunnel at 
a velocity of 100 miles per hour. ' Experimental results are presented 
for angles of attack from -8° to 24° for the swept-wing—fuselage model 
and from -8° to 16° for the unswept-wing—fuselage model. 


The flow characteristics were obtained with a rake of hemispherically 
headed probes utilizing both downwash- and sidewash-angle orifices in con- 
junction with pitot-static orifices to measure dynamic pressure. The 
instrument employed in this investigation is similar to that employed in 
reference 1 and is shown installed on one of the test models in figure 2. 
The flow surveys were made over the right wing with the model inverted to 
minimize support-strut interference and, therefore, represent conditions 
(due to model symmetry) under the left wing of the model. 


Consideration of the angularity rake calibration, data-reduction 
process, method of rake support, possible errors in misalinement, and 
inherent wind-tunnel misalinement angles indicates that the downwash 
data are accurate within approximately $1.0°, the sidewash data are 
accurate within approximately +1.5°, and the dynamic-pressure-ratio data 
are accurate within approximately +0.025. 


THEORETICAL METHODS 


The characteristics of a field of flow can be completely defined 
by the magnitude and direction of the local velocity vectors. It is 
generally convenient to express the direction in terms of the angles e¢ 
in the vertical plane.and ogo in the lateral plane and to express the 
magnitude in terms of local dynamic pressure q;- In order to deter- 
mine the foregoing flow characteristics by use of theory, a knowledge 
is required of the induced velocities contributed by the various 
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surfaces responsible for disturbing the free-stream flow. The discussion 
of the calculative procedure will be restricted in the present section:to 


a brief general description with the specific details and equations enlarged 


upon in appendix A. The principal factors necessary to describe the flow 
sttastlle in #4. z 


characteristics are defined schematically in figure 3. 


In the calculation procedures employed, it was assumed that the flow 
was potential and planar, and, hence, the effects of boundary-layer sep- 
aration and the rolling up and displacement of the trailing-vortex wake 
have been neglected. The effects of the presence of the fuselage have 
also been neglected since the variation of upwash angle induced by the 
circular-cross-section fuselage decays rapidly with lateral distance. 
This variation in upwash angle is presented in figure 4 as a function of 
lateral distance, nondimensionalized with respect to the swept-wing semi- 
span. For the swept-wing configuration, the ratio of fuselage diameter 
to wing span is 0.13. For the lateral locations for which the swept- 


wing calculations have been made, y| 2 = 0.50 and yf 2 = 0.75, the 


fuselage-induced upwash angles are seen from figure 4 to be approximately 
8 percent of wing angle of attack for the inboard location and approxi- 

mately 3 percent for the outboard location. For the midsemispan location 
of the unswept wing, which has a ratio of fuselage diameter to wing span 
of 0.16, the fuselage-induced upwash angle is approximately 10 percent of 


the wing angle of attack. 


The foregoing discussion has considered only the effects of the fuse- 
lage alone. Examination of reference 4 indicates that the mutual- 
interference effects caused by the addition of a wing to the fuselage 
produce only slight changes in the exposed wing-span load distribution. 
Since the calculations of present interest are critically affected by 
lift coefficient and since the comparison of theory with experiment is 
most readily made for comparable lift coefficients, the small changes in 
load distribution indicated by reference } are assumed negligible. For 
regions closer to the fuselage, however, or for larger ratios of fuselage 
diameter to wing span, it is evident from figure 4 that the presence of 
the fuselage should be considered. In this respect, the analyses of ref- 


erences } and 7 may be useful. 


In order to determine the flow characteristics in close proximity 
to the wing, it is necessary to account for both the lift-induced veloc- 
ities and the nonlifting or thickness-induced velocities. The former 
velocities are primarily a function of wing angle of attack and plan- 
form geometric characteristics, whereas the latter velocities are inde- 
pendent of angle of attack and are primarily a function of the local 
airfoil-section thickness distribution, modified by plan-form character- 
istics. Extensive theoretical investigations of the zero-lift velocity 
distributions on the surface of unswept and sweptback wings have been 
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MODELS AND TESTS 


The models about which the flow surveys were made consisted of both 
swept- and unswept-wing—fuselage combinations. Drawings of the wing- 
fuselage combination are presented in figure 1. The wing of the swept- 
wing—fuselage combination had 45° sweep of the quarter-chord line, an 
aspect ratio of 4.0, a taper ratio of 0.3, and NACA 65A006 airfoil sec- 


tions parallel to the plane of symmetry. The wing of the unswept-wing—. 
fuselage combination had 0° sweep of the one~half~chord line, an aspect 
ratio of 3.0, a taper ratio of 0.5, and NACA 65A004 airfoil sections 
parallel to the plane of symmetry. The fuselage consisted of an ogival 
nose section, a cylindrical center section, and a truncated tail cone. 


The fuselage ordinates are presented in table I. 


The tests were made in the Langley 300 MPH 7- by 10-foot tunnel at 
a velocity of 100 miles per hour. ‘Experimental results are presented 
for angles of attack from -8° to 24° for the swept-wing—fuselage model 
and from -8° to 16° for the unswept-wing—fuselage model. 


The flow characteristics were obtained with a rake of hemispherically 
headed probes utilizing both downwash- and sidewash-angle orifices in con- 
junction with pitot-static orifices to measure dynamic pressure. The 
instrument employed in this investigation is similar to that employed in 
reference 1 and is shown installed on one of the test models in figure 2. 
The flow surveys were made over the right wing with the model inverted to 
minimize support-strut interference and, therefore, represent conditions 
(due to model symmetry) under the left wing of the model. 


Consideration of the angularity rake calibration, data-reduction 
process, method of rake support, possible errors in misalinement, and 
inherent wind-tunnel misalinement angles indicates that the downwash 
data are accurate within approximately +1.0°, the sidewash data are 
accurate within approximately $1.5°, and the dynamic-pressure-ratio data 
are accurate within approximately +0.025. 


THEORETICAL METHODS 


The characteristics of a field of flow can be completely defined 
by the magnitude and direction of the local velocity vectors. It is 
generally convenient to express the direction in terms of the angles e¢ 
in the vertical plane. and og in the lateral plane and to express the 
magnitude in terms of local dynamic pressure q,- In order to deter- 
mine the foregoing flow characteristics by use of theory, a knowledge 
is required of the induced velocities contributed by the various 
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surfaces responsible for disturbing the free-stream flow. The discussion 
of the calculative procedure will be restricted in the present section:-to 


a brief general description with the specific details and equations enlarged 


upon in appendix A. The principal factors necessary to describe the flow 
characteristics are defined schematically in figure 3. 


In the calculation procedures employed, it was assumed that the flow 
was potential and planar, and, hence, the effects of boundary-layer sep- 
aration and the rolling up and displacement of the trailing-vortex wake 
have been neglected. The effects of the presence of the fuselage have 
also been neglected since the variation of upwash angle induced by the 
circular-cross~section fuselage decays rapidly with lateral distance. 
This variation in upwash angle is presented in figure 4 as a function of 
lateral distance, nondimensionalized with respect to the swept-wing semi- 
span. For the swept-wing configuration, the ratio of fuselage diameter 
to wing span is 0.13. For the lateral locations for which the swept- 


wing calculations have been made, y/e = 0.50 and yt = 0.75, the 


fuselage-induced upwash angles are seen from figure 4 to be approximately 
8 percent of wing angle of attack for the inboard loration and approxi- 

mately 3 percent for the outboard location. For the midsemispan location 
of the unswept wing, which has a ratio of fuselage diameter to wing span 
of 0.16, the fuselage~induced upwash angle is approximately 10 percent of 


the wing angle of attack. 


The foregoing discussion has considered only the effects of the fuse- 
lage alone. Examination of reference } indicates that the mutual- 
interference effects caused by the addition of a wing to the fuselage 
produce only slight changes in the exposed wing-span load distribution. 
Since the calculations of present interest are critically affected by 
lift coefficient and since the comparison of theory with experiment is 
most readily made for comparable lift coefficients, the small changes in 
load distribution indicated by reference 4} are assumed negligible. For 
regions closer to the fuselage, however, or for larger ratios of fuselage 
diameter to wing span, it is evident from figure 4 that the presence of 
the fuselage should be considered. In this respect, the analyses of ref- 


erences 4 and 7 may be useful. 


In order to determine the flow characteristics in close proximity 
to the wing, it is necessary to account for both the lift-induced veloc- 
ities and the nonlifting or thickness-induced velocities. The former 
velocities are primarily a function of wing angle of attack and plan- 
form geometric characteristics, whereas the latter velocities are inde- 
pendent of angle of attack and are primarily a function of the local 
airfoil-section thickness distribution, modified by plan-form character- 
istics. Extensive theoretical investigations of the zero-lift velocity 
distributions on the surface of unswept and sweptback wings have been 
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reported in references 8 to 11 and indicate that the isobars, that is, 
lines of constant pressure, tend to be parallel to the local lines of 
constant percent thickness for regions not too close to the wing root or 
tip. Reference 9 also shows that the effect of aspect ratio on the back- 
wash velocities is negligible for aspect ratios that are of present inter- 
est (aspect ratios of 4 and 3 for the swept and unswept wings, respec- 
tively). In view of this, and with consideration of the simple sweep 
theory of reference 6, the present paper considers the airfoil sections 
normal to the local lines of constant percent thickness to be two 
dimensional in nature. 


The perturbation velocities of the two-dimensional-airfoil thickness 
distribution may be determined by either conformal transformations as 
reported in references 12 to 14 or by use of the appropriate singularity 
distribution as determined by the methods of reference 5 or 15. ‘The pres- 
ent paper utilized the method of reference 5 in combination with the 
simple sweep theory of reference 6, as described in appendix A, in order 
to account approximately for the effects of either sweep or taper or both. 


In the calculation of the lift-induced velocities, the present pro- 
cedure utilizes, primarily, four horseshoe vortices distributed in the 
chordwise direction at each of 10 spanwise locations, thus making a total 
of 40 horseshoe vortices. The chordwise vortices are assumed to have 
equal circulation strengths but unequal chordwise spacing. The stratagem 
is then to sum the induction effects at points that lie midway between 
any two adjacent chordwise vortices (where possible) for regions near the 
wing chord, and thereby minimize the objectionable singuiarity effects men- 
tioned previously in the "Introduction". This procedure is hereinafter 
referred to as the finite-step method. An illustrative calculation of 
the lift-induced velocities beneath the swept wing is presented in 
table II. 


In calculating the sidewash velocities, the finite-step method 
becomes increasingly inaccurate as the vertical distance from the wing 
chord plane is decreased. Further study of the assumed horseshoe vortex 
system (see appendix A) indicated that the sidewash velocity would approach 
zero as the wing chord plane was approached. This characteristic is not 
consistent with reality in that the lateral gradient in load or vorticity 
implies the existence of sidewash velocities on the wing surface. 


By use of unpublished theoretical studies made by Percy J. Bobbitt 
of the Langley Aeronautical Laboratory (see appendix A), the sidewash 
velocity at the wing chord plane may be estimated and a more realistic 
variation of sidewash velocity with vertical distance effected. 


The velocities induced by a unit horseshoe vortex in the vertical, 
lateral, and longitudinal directions, which are necessary in ‘the present 
methods, were computed by the equations given in reference 16 and are 
presented in tables III, IV, and V for a large range of distances. 
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Comparison of the values predicted by theory with the experimental 
values indicates that the representation of the airfoil-section thick-, 
ness distribution by a two-dimensional singularity distribution (ref. 5) 
modified by simple sweep theory (appendix A) gives excellent qualitative 
agreement for all vertical locations considered. The magnitudes of the 
flow parameters due to thickness are, in general, also well predicted, 
although the downwash angles are underpredicted for the regions immedi- 
ately ahead of the wing chord. 


The flow characteristics at a wing lift coefficient of 0.49 are 
shown in figure 7(b). The chordwise gradients mentioned previously are 
seen to be more severe than for the zero-lift condition (fig. 7(a)). 

For this lift coefficient (0.49) the lift-induced effects, in general, 
completely overshadow the thickness effects and cause large changes in 
the downwash and sidewash angles in addition to reductions in the dynamic- 
pressure ratios. 


Good agreement is in evidence for the downwash angles except for 
the nearest vertical location where the theory overestimates conditions 
immediately ahead of the wing leading edge. This overestimation is pre- 
sumed. to be due to the assumption in the theory of the two-dimensional. 
type of chordwise load distribution that implies full leading-edge suc- 
tion and, hence, unrealistically large induced effects in this vicinity. 


In the case of the sidewash angles (fig. 7(b)), the assumed finite- 
step theory is seen to become increasingly inaccurate as the vertical 
distance from the wing chord plane is decreased. The modified theory 
(see appendix A), which effects a more realistic variation of sidewash 
velocity with vertical distance (particularly near the chord plane), is 
seen generally to agree more closely with the experimental results than 
does the finite-step method. The modified theory was used in the rest 
of the incompressible sidewash calculations presented in this paper. 


The prediction of the dynamic pressures (fig. 7(b)) by use of the 
finite-step method is seen to be good for all chordwise and vertical 
locations presented. : 


Since it has been shown that the decay in the flow distortions can 
be calculated, it would be desirable to consider in more detail the 
predictability of the flow throughout a more complete lift range. A 
comparison of the theoretical and experimental flow fields existing 
15 percent of the local wing chord beneath the midsemispan location of 
the swept wing is presented in figure 8. 


With a change in sign of the flow angles at the most negative lift 
coefficient (Cy, = ~0.53), the conditions existing on the upper or suc- 


tion side of the wing when at positive lift may, because of model sym- 
metry, be examined. The flow parameters indicate the existence of 
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extremely high values of downwash and sidewash angularity as well as 

large dynamic pressures. Examination of the pitching-moment curve pre- 
sented in figure 5 indicates an unstable break at approximately this 

lift coefficient in the positive lift range (Cz, = 0.49), which signifies 

a loss of lift at the wing tip and indicates the existence of nonpotential 
flow. The potential-flow theory utilized cannot then be expected to 


predict the magnitude of the flow parameters for these conditions. 


As the lift coefficient is reduced to Oy, = -0.26, a rather good 


description of the downwash angles is given by use of theory (fig. 8(a)). 
Good agreement is also obtained throughout the positive lift range to 

Cr, = 0.89, which is rather surprising since at this lift coefficient the 
flow on the suction side of the wing is nonpotential. At Cy, = 1.09, the 
theory is seen to overpredict the downwash ahead of the leading edge and 
to underpredict it over the chord proper. This is presumed to be due to 
the rearward movement of the experimental local center of pressure that 
is associated with leading-edge stalling. 


Examination of figures 8(b) and 8(c) indicates that the calculated 
sidewash angles and dynamic pressures are in reasonable agreement over 
the entire lift range with the exception of the extreme cases, (Cy, = -0.53 


and 1.09 where nonpotential conditions exist. 


Tn order to determine the ability of calculations to predict the 
effect of spanwise position on the flow characteristics, a comparison 
with the conditions existing 15 percent of the local wing chord below 
the three-quarter semispan location of the swept wing is presented in 
figure 9. The zero-lift flow angles (fig. 9(a)) and dynamic pressures 
(fig. 9(b)) are well predicted, which indicates that the zero-lift flow 
characteristics are still essentially two dimensional in nature at 


y= -0.75. As the lift coefficient is increased, however, the agree- 


ment between theory and experiment is seen to deteriorate for the down- 
wash angles (fig. 9(a)) in that the theory gives values too high over 
the chord region. This overestimation is presumed to be due to assuming 
a two-dimensional type of chordwise load distribution to exist at this 
spanwise station for (Cy, = 0.23 and to a combination of the aforemen- 
tioned in conjunction with the proximity of the rolled-up tip vortex 
for Cy, = 0.49. In spite of the defects in predicting the downwash 


angles, the sidewash. angles and dynamic pressures are seen to be reasonably 


well predicted. It should be noted that the experimental downwash angles 


are slightly lower at the outboard location (2 = -0.75 in fig. 9(a)) 


than at the inboard location (fg = -0.50 in fig. 8(a)), whereas the 
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sidewash angles are slightly higher. The dynamic pressures appear to be 
relatively unaffected by spanwise station for the two stations presented 
(figs. 8(c) and 9(b)). 


Unswept-Wing Model 


A comparison of the flow characteristics at a distance 15 percent 
of the local wing chord beneath the unswept wing is presented in fig- 
ure 10. The predicted downwash characteristics (fig. 10(a)) are, in 
general, subject to the same discussion and limitations as those for the 
swept wing; the only notable differences were the underprediction of the 
downwash ahead of the leading edge, whereas there was an overprediction 
for the swept wing (fig. 8(a)). The cause of the nonpotential nature of 
the flow above the wing chord plane, as evidenced by the break in the 
pitching-moment curve (fig. 6), is assumed to be due primarily to leading- 
edge separation. 


The comparison between the experimental and theoretical sidewash 
angles below the unswept wing is shown in figure 10(b). As in the case 
of the swept wing, significant chordwise gradients exist under lifting 
conditions. The finite-step theory in which 10 spanwise and 4 chordwise 
horseshoe vortices were utilized is seen to underpredict the sidewash 
angles. Increasing the number of spanwise vortices from 10 to 20 and 
using the estimated surface sidewash velocity (see appendix A) in deter- 
mining the sidewash velocity variation with vertical distance appear 
to provide better agreement with experiment over most of the chord. The 
disagreements existing ahead of the wing-chord leading edge at positive 
lifts are not fully understood, but some of the disagreement may be due 
to support-strut interference effects that have not been assessed. 


The dynamic pressures (fig. 10(c)) appear to be well predicted 
throughout the lift-coefficient range investigated with the exception 
of the largest negative lift coefficient. 


The effects of sweepback cannot be adequately determined throughout 
the lift-coefficient range by comparing the wings of the present investi- 
gation since several geometric differences exist other than the angle of 
sweep. If it is assumed, however, that, for the midsemispan locations, 
the zero-lift flow characteristics are essentially two dimensional, as 
indicated by the ability of two-dimensional theory to. predict the flow 
characteristics, some insight is gained as to the effect of sweep. Com- 
parison of the zero-lift downwash angles and dynamic pressure of the 
swept wing (fig. 8) with the comparable characteristics for the unswept 
wing (fig. 10) indicates that sweep has little effect on these parameters. 
The differences that do exist are felt to be due to the difference in 
thickness ratios. Examination of the sidewash angles (figs. 8(b) 
and 10(b)) indicates that the effect of wing sweep is to induce larger 
sidewash angles, at zero lift, in accordance with simple sweep theory. 


(See appendix A.) 
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Effects of Compressibility 


In the foregoing discussion, the flow-field characteristics were 
for the incompressible case. It would now be desirable to examine briefly 
the effects of compressibility on the flow characteristics. Since no 
experimental data are available at the higher speeds, theoretical com- 
parisons have been made in order to provide at least a qualitative indi- 
cation of the effect of compressibility. 


The calculated compressibility effects, for a subcritical Mach num- 
ber of 0.80, on the flow characteristics at a distance 25 percent of the 
local wing chord beneath the midsemispan location of the swept wing are 
presented in figure 11 for three conditions. The effect of increasing 
the Mach number on the zero-lift flow characteristics is to cause increases 
in both the downwash and sidewash angularities as well as the dynamic- 
pressure ratio, although the basic-flow structure appears to be relatively 
unchanged. In considering Mach number effects for the lifting condition, 
as calculated by the finite-step method, it is convenient to examine the 
effects from two standpoints, namely, the case where a is held constant 
and the case where (y, is held constant. For the constant a case 
(fig. 11), the effect of increasing the Mach number is to cause large 
increases in the positive and negative magnitudes of the downwash angles 
over the complete chordwise range shown and particularly near the leading 
edge. -Large increases in the region of the leading edge are also evident 
in the sidewash angles and large decreases occur in the dynamic pressure 
over the leading-edge portion of the chord; however, the rear 8 percent 
of the chord appears to be relatively unchanged. Some of these effects 
are due to the fact that the wing in compressible flow at constant a 
is generating more lift than the wing in incompressible flow. In order 
to eliminate these additional lift effects, the effects of compressibility 
at constant lift are also presented in figure 11. For this condition, 
the negative and positive magnitudes of the downwash angles are still 
increased over the incompressible conditions. In the case of the side- 
wash angles, however, although the compressible values are slightly higher 
at the leading edge, they are reduced over the chord proper. The com- 
pressible dynamic-pressure ratios still appear to be reduced at the 
leading edge, but to a lesser extent than for the constant a condition, 
and are actually increased beyond the quarter-chord locations. 


CONCLUDING REMARKS 


A theoretical and experimental investigation of the subsonic-flow 
fields beneath swept and unswept wings indicates the existence of signif- 
icant chordwise gradients in the flow characteristics. These gradients 
diminish in severity as the distance from the wing chord plane is increased. 
Increasing the lift coefficient caused large changes in the local downwash 
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and sidewash angles and in the dynamic-pressure ratios. The effect of 
wing sweep at zero lift was to cause increased sidewash angles. 


The theoretical predictions of the flow-field characteristics were 
qualitatively correct in all cases considered, although there were indi- 
cations that the magnitude of the downwash angles tended to be overpre- 
dicted as the tip of the swept wing was approached and that the side- 
wash angles ahead of the unswept wing were underpredicted. 


The effects of compressibility, as calculated by first-order linear 
theory, indicated significant increases in the chordwise variations of 
fiow angles and dynamic-pressure ratios for both the zero-lift and lifting 
eases. The effects of compressibility for the lifting case in which the 
lift coefficient was held constant were less severe than those for the 
constant-angle-of-attack.case. 


Langley Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., April 26, 1956. 
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APPENDIX A 
DETAILED THEORETICAL CONSIDERATIONS 


The purpose of this appendix is to present a more detailed descrip- 
tion of the calculative procedure described briefly in the text. 


The flow is assumed potential and planar, and, hence, the effects 
of boundary-layer separation and the rolling up and displacement of the 
trailing vortex wake are neglected. The effects of the presence of the 
fuselage have been neglected (see fig. 4) for the lateral locations of 


present interest (vf - 0.5 and 0.15). For regions closer to the fuse- 


lage, however, its presence may be considered by methods similar to those 
reported in references 4 and 7. 


A well-established practice in two-dimensional-airfoil theory is 
to consider independently the effects of thickness and the effects of 
angle of attack (ref. 19). The present paper also employs this procedure 
in determining the flow-field characteristics but includes in the non- 
lifting case first-order three-dimensional effects incurred either by 
sweep or taper or both; and in the lifting case, both spanwise and chord- 
wise distributions of vorticity are considered in an approximate manner. 


Nonlifting Case 


In two-dimensional flow, the nonlifting or thickness-induced per- 
turbation velocities are primarily a function of thickness distribution. 
These perturbation velocities, that is, downwash in the vertical direction 
and backwash in the chordwise direction, may be calculated either by con- 
formal mapping techniques, as reported in references 12 to 14, or by use 
of the appropriate singularity (source sink) distribution, as reported in 
references 5 and 15. 


In three-dimensional flow, the problem of determining the perturba- 
tion velocities in the field surrounding the wing becomes considerably 
more complex and requires, in rigorous form, a representation of the wing 
by an infinite number of singularities which must be integrated over’ the 
wing surface (refs. 8 to 11). 


Examination of the extensive theoretical investigations of the zero- 
lift longitudinal or backwash velocity distributions on unswept and swept- 
back wings reported in references 8 to 11 indicated that it is necessary 
to determine only the three-dimensionai effects incurred either by sweep 
or taper or both, since the isobars tend to be parallel to lines of con- 
stant percent thickness (for regions not very close to the wing root or 
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tip) and since the effect of aspect ratio on the local velocities is 
negligible (ref. 9) for the aspect ratios considered in the present 
paper. In view of the foregoing discussion, the following development 
(zero-lift case) will be primarily two dimensional in nature and will 
generally consider swept wings by use of simple sweep theory (ref. 6); 
but the procedure will also be applicable to unswept wings. 


The original contribution of simple sweep theory (ref. 6) was to 
indicate a geometric device by which the critical Mach number of wings 
could be raised. Reference 6 points out that the wing pressure distri- 
bution was chiefly affected by the velocity component normal to the lines 
of constant percent thickness. In determining the zero-lift or thickness- 
induced velocities of a swept wing, it is, therefore, necessary to con- 
sider the thickness distributions of the airfoil sections normal to the 
lines of constant percent thickness. These airfoil sections will herein- 


after be referred to as normal sections in order to differentiate them 
from the streamwise sections. 


The geometric characteristics necessary in the calculation of the 
thickness-induced velocities is shown for the swept wing of the present 
investigation in figure 12. The streamwise chord locations at which the 
flow-field characteristics are desired are indicated bythe data points. 
The normal sections were assumed to be two dimensional and, therefore, 
the perturbation velocities generated by these sections, in conjunction 
with the reduced velocity component V cos A could be calculated by 
either of the two-dimensional-flow techniques mentioned previously (con- 
formal mapping or singularity solution). For the points ahead of the 
wing leading edge, the sweep angles of the normal sections generating 
the perturbation velocities at these points (as indicated by the dashed 


lines in fig. 12) were assumed constant and equal to the sweep angle of 
the leading edge. 


Since the perturbation velocities along and perpendicular to the 
chords of the normal sections (u, and wy, respectively) have been deter- 
mined, it is now necessary to determine the components of these velocities 
relative to the streamwise chord (fig. 12). The downwash velocity w 
remains unchanged since the effects of the increased normal-section thick-~ 
ness ratio relative to the streamwise-section thickness ratio are canceled 
by the reduced normal velocity component. The normal-section backwash 
velocity Uy, must, however, be added to the normal-velocity component 
Vcos A (fig. 12). These vectors are then combined with the parallel- 
velocity component V sin A. This vector addition (fig. 12) determines 
the direction of the resultant-velocity vector Vp relative to the free- 


stream direction. This resultant-velocity direction is seen to be toward 


the plane of symmetry for regions of supervelocity (VR > V) and toward the 
wing tip for regions of subvelocity (VR < V). 


The backwash and sidewash perturbation velocities relative to the 
free-stream direction are (from the vector diagram of fig. 12) 
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U = Un cos A (Al) 
V= Uy sin A (A2) 


and the flow angles in the vertical and lateral directions are, respec« 
tively, 


e = tan-1 W/V tan~+ eee |) ae - (a3) 
ae L+ Un cos A 
V 
Uy sin A 
NN Nore Ste VE (al) 
l+ u 1 Uy cos A 
Vv 


The dynamic-pressure ratios are defined by 


H_ Wu 2 + we + ve (45) 
lo ve 
or, since 
(we + ve) << (V4 u)® 
then 
q 2 Un cos A - 
% (vw? ~ (1, mses (16) 
do ye v | 


In the foregoing development, it was assumed necessary, because of 
wing taper, to determine the thickness distributions of each of the sec- 
tions normal to the lines of constant percent thickness, and then to cal- 
culate the perturbation velocities generated by these sections. It is 
obvious that fulfillment of this assumption would entail a prohibitive 
amount of computational labor. In order to reduce the computations to 
practical proportions, it is necessary to introduce certain simplifying 
assumptions. It was, therefore, assumed that the given tapered swept 
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wing could be replaced by some equivalent infinite-span, swept, untapered 
wing. The effects of wing taper would be retained, however, in using the 
correct local sweep angles in equations (Al) and (A2). 

In order to evaluate the changes in the alrfoil thickness distribu- 
tion incurred by the foregoing assumption, the thickness distributions 
of the normal sections (as indicated by sections 1 to 7 in fig. 12) were 
determined and were found to have maximum thickness ratios of 7.45 to 
7-7 percent. These thickness distributions were then compared with the 
thickness distribution of the streamwise airfoil section which was 
increased. so that its maximum thickness ratio was equivalent to the 
average maximum thickness ratios of the normal sections (7.6 percent). 
This comparison is presented in figure 13. It is evident from this fig- 
ure that wing taper causes some small variations in the thickness dis- 


tributions, particularly over the rear portion of the chord; however, 
when consideration is given to the fact that the maximum surface velocity 
induced on an NACA 65A008 airfoil section is only of the order of 10 per- 
cent greater than the free-stream velocity (for zero lift, see ref. 20), 
it may safely be assumed that these differences in thickness distribu- 


tions, due to wing taper, are negligible. 


Since it has been shown that the given swept wing can be approxi- 
mated by an infinite-span, swept, untapered wing without incurring any 
appreciable differences in the airfoil-section thickness distributions, 
some useful relationships between the assumed infinite-span, swept, unta- 
pered wing and an infinite-span, unswept, untapered wing should be noted. 


Comparison of an infinite~-span, swept, untapered wing with an 
infinite-span, unswept, untapered wing of the same streamwise thickness 
yatio indicates that the normal-section thickness ratio of the swept wing 
is increased by 1/cos A relative to the streamwise section and that 
the normal component of the imposed velocity is decreased by cos A. 


(See the following sketch.) 
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It can, therefore, be reasoned that, since the perturbation velocities 
are linear functions of thickness, for small thickness ratios (as indi- 
cated by an analysis similar to that of ref. 21), the increased thickness 


Th a 
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effects (t = =) are canceled by the reduced velocity V cos A 
we YUP 43) 


perturbation velocities relative to the normal section of the swept wing 
are then approximately equal to the perturbation velocities relative to 
the streamwise section of the unswept, untapered wing; that is, 


Un = Us (AT) 


where ug is the backwash velocity generated by the streamwise thickness 
distribution in two-dimensional flow with a free-stream velocity equal 
to V. 


Equations (Al) and (A2) may now be rewritten as 


u = ug cos A (A8) 


< 
il 


ug sin A (A9) 


and the flow angles given by equations (A3) and (A+) may be rewritten 
as 


eet eet (A10) 
1 Ug cos A 
Vv 
ug sin A 
o = ~ tan7+ = (A11) 
1 Ug cos A 
Vv 
The dynamic-pressure ratio is now 
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$ 
1 The present paper utilized the singularity-distribution method of 
. reference 5 in order to calculate the two-dimensional perturbation veloc- 
ities in the field surrounding the NACA 65A-series airfoils of the swept 
and unswept wings. These velocities were then modified: by the use of equa- 
; tions (A8) and (A9) to account for the three-dimensional-flow effects of 
A either sweep or taper or both. The calculated velocities induced at the 
midsemispan location of the swept wing at zero lift are presented in fig- 
ure 14, and the flow-field parameters determined from equations (A10) 


to (Al2) are presented in figure 7(a) for comparison with experiment. 


Lifting Case 


The general practice of accounting for the wing lift-induced veloc- 
ities, by employing a single lifting line (approximated by a number of 
horseshoe vortices), becomes increasingly inaccurate as the vortices are 
approached. (See ref. 1.) In order to obtain more. realistic values of 
the lift-induced velocities for regions close to the wing, a more detailed 
accounting of the chordwise distribution of vorticity is required. It 
should be noted that, if the actual load distributions are known, they 
would probably greatly enhance the accuracy of the calculations. In the 
absence of these loadings for the wings of the present investigation, the 
spanwise loadings were determined by the method of reference 17 and the 
chordwise load distributions were assumed to be two dimensional in shape 
with the local circulation strength dictated by the span-load distribution. 
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° The shape function of the two-dimensional chordwise vorticity accum- 
: lation Os is given by reference 16 and may be expressed, with a change 
” in variable, as 
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i It was further assumed that this chordwise accumulation could be 
approximated by a finite number of vortices of equal strength since the 

stratagem was to determine where possible, the perturbation velocities, 

1 due to the vortices, at points in the field (in the immediate vicinity 

of the local chord) lying midway between any two adjacent vortex locations, 

thus effecting some cancellation of the objectionable effects of the single 


me lifting line. 
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Integration of equation (Al3) gives the chordwise accumulation of 
vorticity as 


I(x/c)2 


a : = = é) + sin-t \E (Al4+) 


The chordwise limits necessary to insure equal circulation strengths 
(x/ce)3 and (x/c)> must be determined by trial and error. After these 


limits are determined, the centroidal locations of the vortices may be 
found by 


(Al5) 
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A study of the number of two-dimensional-flow vortices needed to 
approximate the airfoil boundary conditions, that is, a= -w/V, in which 
combinations of one, two, four, and eight vortices were considered, indi- 
cated that one and two vortices were insufficient. Utilization of eight 
vortices, of course, was found to give the best approximation of those 
investigated, although this was felt to raise the computations to the 
prohibitive level. Four chordwise vortices were, therefore, chosen as 
the best compromise between required labor and the approximation of the 
boundary conditions. The centroidal locations of these four vortices 
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were found, from equations (Al+) and (A16), to be approximately 
x/e = 0.013, 0.092, 0.272, and 0.621. 

The vortex arrangements thus chosen to represent the wi an for 
consisted of four chordwise horseshoe vortices at each of 10 spanwise 
stations. The vortex arrangement assumed to represent the swept wing 
is presented in figure 15. 


ng plan form 


The equations of the lift-induced perturbation velocities for the 
assumed vortex arrangement may be expressed as 


n=10 m=} 
Ya _ r 
— = —F Al 
Vo oAxVs = mi + " aa 
va_ n=10 m=4 rs (a8) 
Vv UsxVs n=l mt ae 
Wa 1 nel0 m= - 
Vo soAnVs ml mel i nag) 


where Fy, Fy, and Fy, are the geometric functions associated with a 
unit horseshoe vortex. The equations of these functions, as given in 
reference 16, with the appropriate sign changes and nondimensionalized 
with respect to the semiwidth s of the vortex, are presented in appen- 
dix B. The values of these functions over a wide range of distances are 
presented in tables III to V. 


Since 10 spanwise vortices were assumed in the present investigation, 
the semiwidth of each horseshoe vortex is 


b 
— A20 
5 (A20) 


The circulation strength I may also be related to the local section 
lift coefficient by 


(21) 
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Equations (A17) to (A19) may now be expressed as 


Bi. (A22) 


Wa 5 n=10 m=4 cc 
Seg iee E, (a2k) 


The lift-induced velocities were computed for the wing plan forms 
of the present investigation by use of equations (A22) to (A24%) by using 
the span-load distributions presented in figure 16 as determined by the 
method of reference 17. A sample calculation of the lift-induced veloc- 
ities for each unit of lift coefficient for the swept wing is presented 
in table II. The velocities induced at several vertical locations below 
the midsemispan location of the swept wing are presented in figure 17. 


A study of the lift-induced velocities indicated that the downwash 
and backwash velocities calculated by use of equations (A22) and (AQ) 
(fig. 17) had the correct qualitative variation with vertical distance, 
whereas the sidewash velocities did not. Examination of the sidewash 
velocity factor Fy (see eq. (B6)) indicates that when a finite number 
of horseshoe vortices are used the sidewash velocity for small vertical 
distances must approach, at the surface, either zero or become infinite, 
depending on whether the point of interest lies between the trailing 
vortices or directly under a trailing-vortex segment. The points of 
interest in the present calculations were chosen midway between the 
trailing segments of the horseshoe vortices and, hence, approach zero 
as the wing chord plane is approached. In reality, this condition does 
not exist since the lateral gradient in loading or vorticity implies the 
existence of sidewash velocities at the wing surface. Clearly, then, 
sidewash velocities calculated by use of the finite-step method 
(eq. (A23)), where the sidewash velocity is zero at the wing surface, 
would yield much smaller values for points close to the wing (fig. 17) 
than would a method accounting for the finite sidewash at the wing 


surface. 


Unpublished theoretical studies (eqs. (A25) to (Ab2)) made by 
Percy J. Bobbitt of the Langley Laboratory have indicated that a more 
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realistic value of the sidewash velocity variation with vertical distance 
could be obtained by estimating the sidewash velocity at the wing chord 
plane due to the lateral gradient in the velocity potential (referred to 
herein as the chordwise accumulation of vorticity) and then by fairing 

the maximum sidewash velocity in the wing field, as calculated by equa- 
tions (A23) and (B6), to this chord-plane velocity. The sidewash velocity 
at the wing chord plane may be determined from the lateral gradient in 

the chordwise accumulation of vorticity which may be expressed as 


Ve = een) (A25) 


which may be nondimensionalized as 


3 $(x,y) 
VCz, b 
eee . (a26) 
VCcy, y 
b/2 


In the absence of experimental information regarding the chordwise 
accumulation of vorticity ¢g for the wings of the present investigation, 
the two-dimensional vorticity accumulation given by equation (A14) was 
assumed. In order that the total circulation of the system be correct, 
the total chordwise circulation strengths must be corrected to agree with 
the strengths of spanwise vorticity distribution. Thus, equation (Al) 
may be expressed as 


1 abe = ee = (E) + sin! x (427) 
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evaluation of equation (A27) at the trailing edge of the chord 
(x/c = 1.0) gives 


i Ee Gea (a28) 


The three-dimensional vorticity equation given by equation (A421) 
may be nondimensionalized as 


(A29) 


The two-dimensional circulation strength (eq. (A28)) may now be 
corrected to the three-dimensional value (eq. (A29)) by defining a cor- 
rection factor K as the ratio of equation (A29) to (A28). 
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(A30) 
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Multiplying equation (A27) by the correction factor (eq. (A30)) gives 


(A31) 


G(x,y) 1/ Sr 
b rA CL.cay 
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which is the assumed chordwise vorticity accumulation in terms of the 
correct local total circulation strength. 


An approximate expression for the sidewash velocity existing at the 
wing chord plane may now be obtained by substituting equation (A31) into 


equation (A26): 
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Inasmuch as it is difficult to express the geometric characteristics 
of the swept wing in analytic terms amenable for use in equation (A32), 
the required differentiation may best be performed graphically. An illus- 
trated example of this procedure is presented for the swept wing in fig- 
ure 18, and the manner in which the sidewash velocities existing in the 
field are faired to the estimated chord-plane velocity is shown in 
figure 19. 


Further studies of the sidewash-velocity variation with vertical 
distance made by increasing the number of spanwise horseshoe vortices 
also indicated more realistic characteristics except for vertical loca- 
tions very close to the wing chord plane. These characteristics have 
previously been reported in reference 22 for somewhat different circum~ 
stances. The effects of increasing the number of spanwise horseshoe 
vortices on the variation of sidewash velocity with vertical distance 
are shown for the unswept wing in figure 20. 


The flow-field characteristics due to the lift-induced velocities 
may now be determined by 
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Combined Effects 


In order to determine the total flow characteristics, it is nec- 


essary to combine the lifting and nonlifting velocities. The total flow- 
field characteristics may be written as 


w, Wa 
V+ Vo, b 
Bie bans ae (A36) 
Ug cos A Ug 
1+ + —— 
Vv VCq, 


o= - taney (A37) 
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In order to eliminate errors involved in estimating the lift-curve 
slopes of the wings under consideration, the comparisons of theory with 
experiment were made at the same lift coefficient. A comparison of the 
theoretical flow fields with experiment, under lifting conditions, beneath 
the midsemispan location of the sweptback wing as calculated by equa- 
tions (A36) to (A38) is presented in figure 7(b). 
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ffects of Compressibility 


In determining the first-order compressibility effects on the flow- 
field characteristics, the three-dimensional Prandtl-Glauert transforma- 
tion, as given by reference 18, may be used. The general computational 
procedures involved in this transformation have been stated very simply 
by Dr. S. Katzoff of the Langley Laboratory and are presented in the 


subsequent discussion: 


Hb 


The incremental velocities at a point P on the surface of a 
thin body B in compressible flow may be obtained in three steps: 


(1) The x-coordinates of all points of B are increased by 


the factor 1/B, where 8 = \ 1 - Me and where the x-axis is in 
the stream direction. This transformation changes B into a 
stretched body B'. 


(2) The incremental velocities u’, v', and w' in the direc- 
tion of the x-, y-, and z-axes, respectively, at the point P' on 
B' corresponding to the point P on B are calculated as though 
B' were in an incompressible flow having the same free-stream velo- 
city as the original compressible flow. 


(3) The values u, v, and w of the incremental velocities 
at the point P on the original unstretched body B in compres- 
sible flow are then found by the equations 


gadis 

us 32 u (A39) 

v=iy' (40) 
B 

st 5 w (A421) 


It is pertinent to note that the result of step (1), that is, 
stretching the wing chord, causes the transformed wing to have an 
increased angle of sweep, a decreased aspect ratio, a decreased thickness 
ratio, and a decreased angle of attack. The relationship between the geo- 
metric parameters of the given wing in compressible flow and its trans- 
formed equivalent wing in incompressible flow may be expressed as 
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At = raxt(tent (A¥7) 
B 
a’ = Ba (A48) 


The perturbation velocities in the field due to the transformed 
wing in incompressible flow, as indicated by step (2), may now be cal- 
culated by the methods mentioned previously in this appendix. It should 
be noted, however, that, although the chordwise and spanwise locations 
of interest remain unchanged in the transformation, as: indicated by 
equations (A42) and (A45), the vertical locations of interest move closer 
in percent of local chord to the equivalent transformed wing chord plane. 
(See eq. (A43).) 


In accordance with step (3) of Katzoff's general directions, the 
perturbation velocities due to the transformed wing may now be resolved 
into their final form by equations (A39) to (Al1). 


A few specific observations, supplementary. to the foregoing general 


procedure, are appropriate inasmuch as they may somewhat reduce the nec- 
essary computations. 
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Nonlifting case.- If the first step of the transformation, that is, 
stretching. the plan form in the x-direction, which is shown for the swept a 
F wing in figure 21, is assumed to have been completed, it may be observed | 
_ from equation (A¥4) that the thickness ratio is reduced by f. Also, if a 
oe it is noted from equations (A359) to (A4¥1) that the perturbation velocities | 

must be increased by inverse functions of pB, it is apparent that some | 
; beneficial (time saving) cancellation effects might be realized. Care '{ 
ie must be taken, however, that the correct relationship between corresponding : 
vertical locations are used (eq. (A43)). 


In view of the foregoing discussion, it is readily seen that the 
downwash velocity w remains unchanged since the reduced thickness effects 
(eq. (A44)) are canceled by equation (A41). The downwash w at loca- 


tion - ; 2 below the wing in compressible flow is then equal to the 


downwash w ata location -z/¢ below the wing in incompressible flow. 
This simple transformation of vertical locations is possible since the 
downwash velocity at zero lift is independent of the wing sweep angle 
(as shown previously in this appendix). 


i 

\ 

{ 

i 

| In the case of the backwash and sidewash velocities, although some 
| 

{ 

| 

W 


cancellation of the thickness effects are realized, a simple transforma- 
tion of vertical distances is not immediately possible since these velo- 
cities are also a function of the transformed wing Sweep angle (eqs. (A8), 
(a9), and (A47)) Some saving is possible, however, by considering equa- 
. tions (a8), (39); (A39), (A40), and (A#7), and noting by use of equation 


(A44) that ug' = Bug, from which the following may be deduced: 


| V= ug sin see (A49) 


i Ug cos A cos A' 


A wae (A50) 
hoo B cos A 


\ where again the corresponding vertical locations in compressible and 
: incompressible flow (as given by eq. (A43)) must be observed. 


With the perturbation velocities now determined, the flow-field 
4 characteristics in compressible flow, for subcritical Mach numbers, for 
. nonlifting conditions may be found by equations (Al0) to (A112). 


The calculated first-order zero-lift compressibility effects, for 
a subcritical Mach number of 0.8, on the flow-field characteristics 
beneath the midsemispan location of the swept wing are presented in 
figure 11. 
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Lifting case.- In calculating the effects of compressibility on the 
lift-induced perturbation velocities, it is necessary to follow only the 
general outlined procedure. The perturbation velocities at corresponding 
vertical locations (given by eq. (A43)) may then be expressed, by use of 


spe e A LAAN we ‘ 
equations (A@2) to (AG4) and (A39} to (All), as 


yer 


Ug 1 U,' 

—_F- = A51)- 

VCy, 8° Vc, (451) 
' 

vc; «6B vC;,' 

fog (53) 

vcy, 8B VCy," 


If comparing the effects of compressibility on the flow-field char- 
acteristics on a constant a basis is desirable and the calculations 
‘are performed on the basis of unit lift coefficient, as it is generally 
convenient to do, some care must be exercised in the lift-coefficient 


reduction in order to obtain the proper a. 


Since 
On! = (Cry) 'a! (A5}) 
then substituting equation (A48) into equation (A54) gives 
Oy,’ = (Cr, )'B« (A55) 
where (CL.)' is the lift-curve slope of the equivalent transformed 


wing and is not to be confused with the true compressible lift-curve 
slope. 
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The equations for the perturbation velocities (A51) to (A535) for a 
constant a comparison may now be expressed by 


= ile Ua" t 
—_z = : Cc a (A56) 
(= a=Constant B VC, ( Te.) 

*s) Vi : 

= = ——(C,,,)'o (A57) 
( V /g=Constant VCy, ( a) 

*s) We 

= eee ae be (A58) 
( a=Constant VCy, ("tx) 


The calculated compressibility effects, at constant a, on the flow- 
field characteristics beneath the midsemispan location of the swept wing 
calculated by the aforementioned equations and combined with the zero- 
lift perturbation effects are presented in figure 11. 


If it is desired to determine the calculated effects of compres- 
sibility on the flow-field characteristics on the basis of constant lift 
coefficient, it is necessary to decrease only the lift-induced perturba- 
tion velocities at constant a, as given by equations (A56) to (A58), by 
the ratio of the incompressible lift-curve slope to the true compressible 
lift-curve slope. 


The compressible lift-curve slope of the swept wing used in the 
present paper was determined from the equation 


(A59) 


This expression, which was developed by Edward C. Polhamus of the 
Langley Laboratory in 1949, is an improved version, with regard to low 
aspect ratios and compressibility effects, of that presented in refer. 
ence 23. Another, but somewhat more complicated, form of this equation 
has been independently developed in reference 24. With regard to the 
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use of the sweep of the half-chord line in equation (A59) , a recent 
unpublished analysis by Polhamus indicates that there is little effect 
of taper ratio for wings having the same half-chord-line sweep angles. 


The calculated compressibility effects, at constant lift, on the 


flow-field characteristics beneath the micseu epee location of the swept 


wing are presented in figure 11. 
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The positive directions of distances and velocities used in deter- 
mining the induction characteristics of a unit horseshoe vortex are 
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DOWNWASH 
i 

, 

: defined in the following sketch: 
i: 

' 

\ 

i 
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Downwash Equation 


x/s 


The downwash velocity induced at a point in space is given by the 


following equation: 
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Some identities, due to the symmetry of the aforementioned equations, 
which increase the useful range of table III are given by 


ry(S, M, M2) = vy(is, - 4, £2) 
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Sidewash Equation 


The sidewash velocity induced at a point in space is given by the 
following equation: 


oes Fy (B5) 


where 
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Some identities, due to the symmetry of the aforementioned equations, 
which increase the useful range of table IV are given by 


(B7) 


and. 
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Backwash Equation 


The backwash velocity induced at a point in space is given by the 
following equation: 
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Some identities, due to the symmetry of the aforementioned equa- 
s tions, which increase the useful range of table V are giveri by 
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TABLE IT.- SAMPLE CALCULATION OF LIFT-INDUCED VELOCITIES 
BENEATH THE SWEPT-WING MODEL BY USE OF 
EQUATIONS (A22) TO (AQ) 


{fo ska Pa a 
I¥/3- ~0.53 B= 0.45; 2 = “0.20 


-0.45168 
~ -38099 


~.2HBOL 
110925 


FUN | FUNDH | FUND 


~0.02009 
~ 02003 
- «01969 
~-01733 


-0.01851 
ee 
-.0175' 

~-01468 


-0.02437 
- 02377 
02186 
~ 01824 


>@® = -0.7172 


Whe vertical distance z/c = -0.10 is identical with Az/s = -0.5 and is eonstant for this table. 
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. , Figure 3.- Sketch showing coordinate system and positive directions of 
velocities and angles. 
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Figure 6.- Lift, drag, and pitching-moment characteristics of the unswept- 


wing—fuselage configuration. 
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(a) Downwash angles. 


13 


Figure 8.- Flow characteristics at the midsemispan location of the swept 


wing for various lift coefficients. z/c = -0.15. 
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(b) Sidewash angles. 


Figure 8.- Continued. 
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(c) Dynamic-pressure ratios. 


Figure 8.- Concluded. 
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(a) Downwash and sidewash angles. 


z/e = -0.15. 


Figure 9.- Flow characteristics at the three-quarter semispan location 
of the swept wing for various lift coefficients. 
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(b) Dynamic-pressure ratios. 


Figure 9.- Concluded. 
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(a) Downwash angles. 
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swept 


pan location of the un 
z/ce = -0.15. 


wing for various lift coefficients. 


Figure 10.- Flow characteristics at the midsemis 
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(c) Dynamic-pressure ratios. 


Figure 10.- Concluded. 


Sa Rl al Nguyen 
. 


Q 
PIS TT 


HEHE HGz = 5 4 | 


PP Eee A TT 
PT MELT TTT IY TT 
PIQELL ET TT er 


nn re eee ae a, 


—ah 


cap enmenggy ene er a ee eee en ee ee en 


1H 
NACA TN 3738 81 
M 
O 
i= ip8 
-4 
| 


€,de9 ¢. Ly 


~4 


Figure 11.- Calculated effects of Mach number on flow characteristics 
beneath the midsemispan location of the swept wing. z /c = -0.25. 
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Figure 12.- Geometric characteristics of wing used in simple sweep theory. 
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Figure 15.- Vortex arrangement assumed to approximate swept-wing lift 
characteristics. 
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(a) Swept wing. 
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Figure 17.- Calculated additional velocities at the midsemispan location 
of the swept wing for unit 1ift coefficient. 
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Figure 18.- Schematic illustration of graphical differentiation to deter- 
mine sidewash velocity on chord plane of swept wing. 
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Figure 19.- Variation of sidewash velocity with vertical distance below 
swept wing. x/c = 0.20. 


90 NACA TN 3738 


——— Vortices at /O spanwise 
locations. Equations (A 23) 
and (B6) 


—— Vortices at 20 spanwise 
locations. Equations (A 23) 


and (B6) 
aN 

Ser 
Vo MESSE RE 

MO, GT estimated chord-plane - 
_2 a jelocit WAS /Ag3e, hee ate ct 
Q «10 -20 -30 -40 +-50 -60 

“te 


Figure 20.- Effect of number of spanwise horseshoe vortices on sidewash 
velocity variation with vertical distance beneath the unswept wing. 
x/c = 0.10. 
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